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Abstract 

We propose a solvable model of a one-dimensional harmonic oscillator quantum gas 
of two sorts of particles, fermions or bosons, which allows to describe the formation 
of pairs due to a suitable pair interaction. These pairs we call "pseudo-bosons" 
since the system can be approximated by an ideal bose gas for low temperatures. 
We illustrate this fact by considering the specific heat and the entropy function for 
N = 8 pairs. The model can also be evaluated in the thermodynamic limit if the 
harmonic oscillator potential is suitable scaled. 
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1 Introduction 



It seems to be generally taken for granted that a compound of an even num- 
ber of fermions behaves like a boson since any state of 2nN fermions has the 
same permutational symmetry as a state of iV bosons. Accordingly, the re- 
cent experiments with trapped ultracold gases [1] [2] [3] [4] are viewed as an 
achievement of Bose-Einstein condensation [5]. However, a moment's thinking 
reveals that this cannot be true in an exact sense but at most approximately. 
Consider an ideal gas of N "bosons" , each "boson" consisting of 2n fermions 
bound together by some attractive force. In the high temperature limit the 
bounds between the fermions will be broken and the system behaves like an 
ideal fermi gas at high temperatures. Its energy will thus be E ~ 2nN fpksT ', 
where fp is number of degrees of freedom for a single fermion, e. g. fp = 6 
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for fermions sitting in a 3-dimensional harmonic oscillator potential which we 
assume now for sake of simplicity. This would be consistent with the corre- 
sponding formula for a bose gas of iV particles if the number of degrees of 
freedom for a single boson is assumed to be f B = 2nf F . For intermediate tem- 
peratures T which are large with respect to the spacing of 1-particle energies, 
say ksT 3> hw, but small with respect to the binding energy of the fermionic 
compounds, we would expect a "preliminary" classical limit of the energy, say 
E ~ 6NksT, different from E ~ 12nNksT. Hence the energy function E(T) 
of the considered gas cannot coincide with that of an ideal bose gas at medium 
and high temperature. 

More interesting with regard to Bose-Einstein condensation is the low tem- 
perature limit. N non-interacting bosons would occupy the same 1-particle 
ground state for T — > 0, which is impossible for 2nN fermions according 
to the Pauli principle. Hence we expect a non-bose behaviour for "bosons" 
also for low temperatures. In view of these problems we have scruples about 
identifying compounds of 2n fermions with bosons and will rather call them 
"pseudo-bosons" . 

To be sure, there will be a close resemblance between bosons and pseudo- 
bosons and the thermodynamic functions of ideal gases of each sort should be 
approximately the same for some range of temperature, not too low and not 
too high. To make this statement more precise is a challenge to the foundations 
of quantum statistical mechanics, in particular, to demonstrate the possibil- 
ity of pseudo-Bose-Einstein condensation. We do not know of any attempt to 
tackle this problem, with the exception of an estimate of the critical density 
of a pseudo-Bose-Einstein condensate where fermionic effects would occur [6]. 
As a first step to treating the problem of pseudo-bosons we seek for simple 
models describing the formation of pairs of fermions, being the simplest case 
of pseudo-bosons, which are solvable in the sense that the partition function 
of the system can be explicitely calculated. For this purpose the model need 
not be a realistic description of, say, rubidium atoms in a trap. It rather serves 
as an example to study the problem of pseudo-bosons in principle. 
In section 6 (summary and outlook) we will come back to the problem of 
pseudo-bosons. Of course, the model we are going to present can be consid- 
ered in its own sake without reference to the specific problem and could be 
compared with other 1-dimensional solvable models of quantum systems, see 
[7]. It turns out that the case of pairs of interacting bosons is completely 
analogous and can be treated simultaneously, although the discussion of the 
results will be confined to the fermi case. The model will be defined in section 
2 and consists of two sorts of iV fermions (resp. bosons) with equally spaced 
non-degenerate 1-particle energy levels. In this aspect it remotely resembles 
the Luttinger model [8], [9], which is, however, unitarily and hence thermody- 
namically equivalent to a system of non-interacting particles. In our model the 
interaction is introduced in such a way that pairs of fermions of different sort 
with the same quantum number gain a binding energy V. Hence the diago- 
nalization of the total Hamiltonian is trivial. The calculation of the partition 
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function in section 3, however, requires some combinatorics. For N < 8 we 
calculate the specific heat, using computer- algebraic means, and discuss its 
properties (section 4). For larger iV one has to resort to the thermodynamic 
limit which is derived in section 5. The convergence to this limit is shown to 
be rather rapid for T > 0. 



2 Definition of the model 



We will consider two sorts of particles with Hilbert spaces and for 
1-particle states. Let ((f)^ eN and (0^) gN be orthonormal bases in 

and H.( 2 \ Define the 1-particle Hamiltonians by 

ff<»fy« = nfr^g\ i = l,2, (1) 

which holds for example if all particles are assumed to sit in a 1-dimensional 
harmonic oscillator (HO) potential and the zero point energy is neglected. 
Usually we will set hio — 1 in what follows, but later we consider a scaling law 
uj = u(N) in the context of the thermodynamic limit iV — > oo. 

The interaction between particles of different sorts will have the simple form 

V^^®^ = -V^®^8 nm , (2) 

where V > is a real parameter. Thus the systems gains the binding energy V 
if two particles of different sorts have the same quantum numbers and form a 
compound, say, a "molecule" . The molecules can be viewed as pseudo-particles 
in a HO potential with double spacing 2fkv of the energy values. 

We consider N fermions (resp. bosons) of each sort, but particles of different 
sorts are still assumed to be distinguishable. The total Hilbert space will thus 
be ± 

n ± = ) ® (®n (2) ) , (3) 

where the ± sign stands for the completely symmetric (+), resp. antisymmetric 
(— ) subspace. The model is defined by the Hamiltonian 

N N N 

tf = E^ 1} + E#f + E v^, (4) 

where the subscripts //, v indicate the factor of the tensor product where the 
1- and 2-particle Hamiltonians of the form (1) and (2) act upon. Clearly, (4) 
is invariant under permutations n G Sn x Sn and will be understood to be 
restricted to the subspace H ± . 
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In the fermionic case the eigenstates of H can be chosen as the products of 
Slater determinants 

N N 

S nm = A ® A < , (5) 

i=l j=l 

where all quantum numbers of the same sort are strictly ordered, ri\ < n 2 < 
. . . < tin and mi < m 2 < . . . < mjy. We will write the latter condition in the 
compact form nm e I - . Analogously one obtains the eigenfunctions in the 
bosonic case, but here we have n\ < n 2 < . . . < un and mi < m 2 < . . . < mjv 
or, equivalently, nm e J + . The corresponding energy eigenvalues are 

N N N 

E nm = J2 n i + J2 m j~ V J2 ^rurrij ■ (6) 
i=l j=l i,j=l 



3 The partition function 

The above defined model is admittedly a rather crude description of a real 
quantum gas, but it allows an explicit calculation of its partition function 
without any further approximation. As usual, the partition function is defined 
by 

Zn(P)= E exp(-/?£ nm ). (7) 

nme/± 

In view of (6) we will split the series (7) into a finite number of subseries 
according to the number of coincidences between the two sequences of quan- 
tum numbers n and m. More precisely, we will split the total index set ^ 
or, equivalently, the set of eigenstates of H into disjoint subsets according to 
their "order type". The order type of nm comprises the linear ordering of 
the union of both sequences n and m, including the coincidences of quantum 
numbers. For N = 2 and the fermionic case we have 10 different order types, 
for example (rii < n 2 < m x < m 2 ) or {n\ = m\ < m 2 < n 2 ). 

In general, every order type can be uniquely characterized by a string s com- 
posed of symbols Sj G {0, 1, 2}. "sj = 1 " means "a quantum number of sort 1", 
analogously for "s, = 2 " , but "sj = " means that "two quantum numbers of 
different sort coincide". In the above example, {n-y < n 2 < mi < m 2 ) <-> (1122) 
and (rii = m x < m 2 < n 2 ) <-> (021). 

Let iVj(s) denote the number of occurences of the symbol "i" in the string s 
and \s\ the length of s. We define ^4^ as the set of all finite strings s over the 
alphabeth {0, 1, 2} such that i\q(s) = N 2 (s) = N - iV" (S). Then the strings s 
in ^.at are in 1 : 1 correspondence with the order types, denoted by T ± (s), of 
the subsequences of (7). Hence we may write 
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Zn(P)=T, E ffi P(-^n) (8) 
seA N mneT±(s) 

= E ^ E E exp (-0 f £ * + E ^) ) • (9) 

N =0 seA N , NQ mneT±(s) \ \ i j J J 

In the last equation we have used the definition 

An,n = {se A N \N (s) = N } (10) 

and that all strings in An,n lead to the same interaction energy —N V. 

Note that the series DmneT±(s) is on ly a 2iV — iVo-fold series because of the N 
coincidences between the quantum numbers. In the exponent, however, each 
term in J2i n i + J2j m j has to be counted once, irrespective of coincidences. 
Hence we may introduce 2N — N new summation indices Oj if we account 
for the double occurences in the exponent by factors gi. The condition a\ < 
ai < ... < a\ s \ will be written in the compact form a e J~ , analogously 
oi < a 2 < . . . < a\ s \ «aG J + . 

This yields 

exp -/? + J2 m i ))= E ex P -/^E^( s ) a i )> ( n ) 

mn€T±(s) \ \ i J / / aeJ± V ' / 

where we have used |s| = 2iV — N and the definition 

»(.) = (i s, = ,; or2 - (12) 

2 Si = 



The remaining evaluation of (11) is straight forward. We introduce as new 
summation variables the differences 



di = a\ (13) 
d n = a n - a„_i, n > 1 (14) 

and use the partial sums 

M 

G n (s) = ]>>(*) (15) 

1=71 

since 

£ G n <i n = £ G n (a n - a„_i) = XX G ™ ~ C n+ i)a n = £ £f n a n . (16) 

ra n n n 
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Further let 



Then we obtain 



£ ±_ r for bosons(+) 
for fermions(— ) 



oo oo 



Eexp -/? =E E ••• E exp -/? E^H 

aeJ± V \ i // d 1 =0d 2 =e ± (i| s |= £ ± V V i // 

oo oo oo 

= ^ e -/3Gidi ^ e -/3G 2 d 2 _ _ _ e _/3G i{if8) 

(ii=0 (i2=£ ± d| B | =e ± 

1 |s| - £ ±/3G„ 



1 _ e -HG_ 
1 C n=2 



The result can be further simplified by utilizing the fact that (19) depends only 
on the G n , i. e. on the positions of the symbol "0" in the string s. There are 
exactly ( 2 ^T^) strings in An,n with the same positions of "0" and hence 

the same values (19). We may thus insert the factor ( 2 ^Tw^) and rather 
sum over the set Bn,n of all strings s over the alphabeth {0, 1} satisfying 
iVo = iVo(s) = N — Ni(s)/2. gi(s) and Gi(s) are defined analogously as before. 
Then the final result reads 

^(«=E^ w 2( ;_* ) £ T^n^^y. (20) 

^0=0 \ iV iV ° / seBjv,iv n=2 1 e 



4 Specific heat and entropy 



As usual, we define 

C±(/3)=/3 2 ^lnZ±(/3). (21) 

The explicit form would be too complex to be reproduced here. However, using 
a computer algebra software the specific heat function can be calculated and 
plotted for small N, say N < 8, without problems. A typical plot of as a 
function of i? = ^£ = is shown in figure 1. 

For i? — > the specific heat vanishes exponentially as exp ( — j^), which is 
typical for finite quantum systems. For i? a first plateau of the height 

Cg = 8 is reached which is a "preliminary" classical limit of a system of 8 
molecules moving in a HO potential, in accordance with the law of Dulong- 
Petit. In the region < d < $i the molecules are practically never broken into 
parts and behave as pseudo-particles. This is confirmed by a comparison with 



6 




Fig. 1. The specific heat C% as a function of i? = with V = 100. This can be 
compared with the specific heat of system of 8 bosons in a 2fkv-HO potential or of 
16 fermions in a hu-KO potential. 

the specific heat function of 8 bosons (or fermions) in a 2huj-HO potential, 
see figure 1. (For the identity = in 1-dimensional HO potentials see 
[10].) For i?>i?! the specific heat increases above its preliminary classical 
level which indicates a beginning dissociation of the molecules. For d ~ $3 the 
final classical limit C$ = 16 of 16 particles in a HO potential is reached and all 
molecules are dissociated. But before, at 1? ~ i? 2 , a marked peak of the specific 
heat occurs which means that at this temperature the maximal portion of any 
energy supply would be used for the dissociation of the molecules instead of 
heating up the system. This peak at 1? ~ $2 might indicate a diffuse phase 
transition, however, the question whether a real phase transition occurs can 
only be addressed by considering the thermodynamic limit N — ■> 00 in the 
next section and will be answered negatively there. 

The entropy function S(T) interpolates between the entropy of 8 ideal bosons 
for low temperature and of 16 fermions of two different sorts for high temper- 
ature, see figure 2. 
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Fig. 2. The entropy S as a function of $ = -jjy,V = 100, for 8 fermion pairs. 
Compare the corresponding plots of the entropy of an ideal gas of 16 fermions of 
two different sorts for high temperature and of 8 bosons for low temperature. 

5 The thermodynamic limit 

It turns out that the explicit solution (20) can also be evaluated in the limit 
iV — > oo. But we have to perform a simultaneous scaling of the energy level 
spacing ftw. This is analogous to the usual practice of letting the volume of a 
system go to infinity in such a way that the particle density remains finite. 

To find the correct scaling law in our case, we consider N fermions occupying 
the first N eigenstates of a HO potential (the "Fermi sea"). Their 1-particle 
density assumes values which diverge as \H\f within an interval of length 
~ \fN in accordance with the total integral of the 1-particle density being N. 
The typical unit of length of the HO is xq = \f^- Hence a scaling uj — > ^ 

stretches the system by a factor \f~N and reduces the 1-particle density to 
finite values for N — > oo. For a more detailed account see [11]. In the case of 
a boson gas we cannot argue in the same way as above, but we will choose 
the same scaling. In this case the existence of a thermodaynamic limit, which 
is shown below, is the only justification of choosing the scale transformation 

U ~* N ■ 
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In practice the considered scaling means that the variable (3 on the rhs of (20) 
has to be replaced by except for the factor e^ N ° v which comes from the 
interaction, not from the HO potential. Consequently, the corresponding terms 
can be approximated by their high temperature limit, i. e. by their first order 
Taylor expansion with respect to (3/N. In this limit the difference between the 
statistics of fermions and bosons disappears. 



Hence 

e - £ ±£ G „( s ) N 



1 _ e -£G„w (3G n {s) 

and we obtain 



(3 fixed, N -> oo (22) 



=£--'(T:3(f)"„E,ij,i « 



In the next step we calculate the combinatorical factor: 



M x x 
J£ w „ Si ^ = ^N \(2N-2N )\ ■ (25) 

In order to prove (25) we first consider an example. Let N = 6, Nq = 3, N\ = 6 
and consider the string s = (101100111) G Bn,n - We then have (g n ) = 
(121122111) and (G n ) = (12,11,9,8,7,5,3,2,1). Hence R n ± = ±2|i. Gen- 
eralizing this example, we see that Jin ^~ can always be written as a fraction 
with denominator (27V)! and the nominator being a product of N integers n u 
with 1 < n v < 2N such that their distance is \n v — n M | > 2 for v ^ fi. Con- 
versely, any N integers with these properties occur in the nominator of the 
terms in the sum Y,seB N Nq Hn=i We set 2N = M and generalize to arbi- 
trary natural numbers M. Let F(M, N ) be defined as the sum of all products 
of iVo-tuples (rii, . . . , njv ) with n\ < n 2 < . . . < njv < M and \n u — n u+ i\ > 2. 

If we enlarge M by 1 additionally those iV -tuples satisfy the above condition 
with ri/Vo = M. Hence the following recursion relation holds for F: 

F(M + 1, iV ) = F(M, N ) + MF(M - 1, N - 1) . (26) 

Moreover, 

F(M, N ) = for M < 2N (27) 

and 

F(M, 0) = 1 for all M . (28) 
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The recursion relation (26) and the initial conditions (27), (28) are satisfied 
by the following solution 

F W-™ = 2*NAM-2N,)V (29) 

as can be easily confirmed. Setting M = 2TV and dividing by (27V)! completes 
the proof of (25). 

Inserting (25) into (24) yields 

N /N\ 2N ~ No 1 

(?) w-*,),. - < 3 °> 

Next we approximate the factorials in (30) by Stirling expressions, e. g. 7V ! ~ 

NQ°e~ N ° , set 7V = TVrr and replace the sum J2n =o • • • by an integral Jq N dx 

The logarithm of the integrand I(x) reads 

In I(x) = TV {2 - x + V(3x - x In 2 + 2(x - 1) ln(l - x) - x In x - 2 In (3 + x In (3) 

(31) 

The term in brackets has as a maximum at 

x = -^- yfi (l + Pe vfi - yjl + 2(3e^ , (32) 

hence the integrand will be sharply peaked at this maximum if TV — > oo. Thus 
we can evaluate the integral for TV — > oo using the saddle point method. 

In the last step we remove all factors in the partition function approximation 
which are not of the form ( N , ( being independent of TV. These factors would 
only give finite TV corrections in the thermodynamic functions and have to be 
neglected in the thermodynamic limit. The final result is 



ZM = 2- 2 ^exp ||e-^ [f^ - 1 + y/l + 2e^/?) } 



-2N 



1 + ^1 + 2e v ^ 



(33) 

The resulting specific heat per molecule will be independent of TV: 



<P)=TF^Zoo(P) ( 34 ) 



Nd(3 2 

2(2(1 + yjl + 2e v ?(3) + e V(3 (3{h + 2V (3 + V 2 (3 2 + y/l + 2e v P(3)) 

(35) 



1 + 2e v ^(l + J I + 2e y Pf3) 
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Fig. 3. The specific heat per molecule C^/N as a function of i? = ^7 with V = 100 
for N = 2, . . . , 8 which rapidly converge to the thermodynamic limit function c(f3), 
except for low temperatures. The tilde in C^/N indicates the scaling to — > ^. 



The rapid convergence of the finite iV specific heat per molecule, scaled with 
uj — > jf, to the limit (35) is illustrated in figure 3. 

The limits iirng_ >0 C (P) — 2 and Hindoo c{(5) = 1 are in accordance with our 
discussion in section 4. 

We define the "degree of dissociation" v(f3) as the thermal expectation value 
of 1 — = 1 — x, or, equivalently, as 



v ^ = 1 -^pW lnZoom (36) 

and obtain the simple result 

u(P) = 2 . (37) 

1 + ^1 + 2^ 

A plot of the degree of dissociation as a function of 1?, see figure 4 confirms 
our previous interpretation of the specific heat function. 

Since (33) is a smooth function of (3 there is no phase transition in our model. 
Instead of a critical temperature at which the molecule phase is changing 
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Fig. 4. The degree of dissociation v as a function of i? = -jjy,V = 100, in the 
thermodynamic limit. 

into the fermi gas phase we have an extended temperature domain where the 
transition occurs, even in the thermodynamic limit. The reason of this may 
may be that our model is only 1-dimensional and, moreover, is not suited to 
describe long range order. 



6 Summary and outlook 

We have studied a two-species quantum gas of 2N particles, fermions or 
bosons, confined in a 1-dimensional HO-potential with a simplified attractive 
2-body interaction between particles of different sorts. The canonical partition 
function and hence the relevant thermodynamical functions of this system can 
be written in terms of finite sums and products and explicitely calculated by 
using computer-algebraic means for N < 8. The specific heat function shows a 
"preliminary" classical saturation value of c ~ Nks for medium temperature, 
a peak indicating the increasing dissolution of the pairs and a high temperature 
limit of c ~ 2Nk B . Using the scaling uj — > ^ we calculated the thermodynamic 
limit, i. e. the limit of lnZj^ and related functions for N — > oo. In this limit 
the difference between bosons and fermions vanishes. The peak of the specific 
heat function remains smooth and no phase transition occurs. 
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We motivated the construction of our model by considering the problem of 
pseudo-bosons, that is, the problem, in which sense pairs or, more generally, 
even numbers of fermions could be considered as bosons. What can be learnt 
from our model concerning this question? As expected, the behavior of the 
specific heat function mentioned above deviates from that of an iV-particle 
bose gas, see figure 1, similarly for the entropy function, see figure 2. 

However, for low temperatures the difference between bosons and pseudo- 
bosons is not so marked as we would expect in the general case. This is partly 
due to the circumstance that the specific heat functions of ideal fermi and 
bose gases in 1-dimensional HO potentials are the same, see [10]. Only the 
ground state energies differ by a constant. Also the entropy function S(T) of 
the pseudo-bose gas is very well approximated by the corresponding bose gas 
function at low temperatures. We expect that the differences would become 
more marked in the 3-dimensional HO potential case, where Bose-Einstein con- 
densation occurs. Unfortunately, it is not clear how to generalize the method 
used in this article to that case, and the problem of pseudo-Bose-Einstein 
condensation remains unsolved for the present. 
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